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Abstract 

An investigation is performed to evaluate the performance of a computational fluid dynamics (CFD) 
tool for the prediction of the reacting flow in a liquid-fueled combustor that uses water injection for 
control of pollutant emissions. The experiment consists of a multisector, liquid-fueled combustor rig 
operated at different inlet pressures and temperatures, and over a range of fuel/air and water/fuel ratios. 
Fuel can be injected directly into the main combustion airstream and into the cavities. Test rig 
performance is characterized by combustor exit quantities such as temperature and emissions 
measurements using rakes and overall pressure drop from upstream plenum to combustor exit. 
Visualization of the flame is performed using gray scale and color still photographs and high-frame -rate 
videos. CFD simulations are performed utilizing a methodology that includes computer-aided design 
(CAD) solid modeling of the geometry, parallel processing over networked computers, and graphical and 
quantitative post-processing. Physical models include liquid fuel droplet dynamics and evaporation, with 
combustion modeled using a hybrid finite-rate chemistry model developed for Jet-A fuel. CFD and 
experimental results are compared for cases with cavity-only fueling, while numerical studies of cavity 
and main fueling was also performed. Predicted and measured trends in combustor exit temperature, CO 
and NO x are in general agreement at the different water/fuel loading rates, although quantitative 
differences exist between the predictions and measurements. 


Introduction 

Natural-gas-fired, water-injected industrial gas turbine combustors have been operational for many 
years (e.g., GE-LM2500), with more recent efforts focused on inlet fogging. Water injection into the 
combustion airstream provides enhanced efficiency with decreased emissions at lower costs. For example, 
dropping the compressor inlet air temperature 11 °C (20 °F) of a 100-MW facility, by fogging the inlet, 
produces nearly a 1 0-percent increase in power, an overall turbine efficiency increase of 3 percent with 
30 percent less fuel, 20 percent less CCL and 10 percent less total NO x (Peterson, 2003). Compare capital 
costs of $20/kW for inlet fogging that can be installed in 3 months to new plant costs of $750/kW that 
require 3 years to build (Peterson, 2003). Further, the GE-7EA turbine as normally steam injected and 
evaporatively cooled producing about 140 ppm NO x , with inlet fogging will experience an additional 
1 8-percent reduction in NO x (Mee, 1999). 
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These benefits may not be entirely realized in aircraft gas turbomachines (or liquid-fueled power 
generation turbines), yet the B-47 with six J47-GE-25A turbojets at 6000 lb of thrust (dry), or 7200 lb of 
thrust (using water injection for takeoff) provided a 20-percent increase in thrust necessary to get off the 
runway. Still today, B52-B, powered by eight PW J— 57— 19 turbojet engines, produce 12 000 lb of thrust 
with water injection, and the water-injected 747-200 is an old yet well proven system (Daggett, 2004). 

In a preliminary emissions study of the impact of water-injected systems on design, airplane 
performance, maintenance, and cost implications for takeoff and climb-out only, (Daggett, 2004) found 
water misting the low-pressure compressor may be preferable to either direct combustor injection or using 
the high-pressure compressor. If the water misting rate were increased from current industrial gas turbine 
rates of 0.83 percent water-to-air ratio to near 2.2 percent, NO x emissions could be reduced some 
47 percent from nonwater misted engines without specific fuel consumption (SFC) penalties and 
estimated reductions up to 242 °C (436 R) in turbine inlet temperature, which could lead to increased hot 
section life. A 0.61-db noise reduction was also calculated. Daggett estimates a NO x emissions reduction 
cost/benefit ratio of about $ 1,663/ton. 

More in-depth studies are required to substantiate these benefits. In one such study, Eames (1988) 
determined that a water- and methanol-injected engine for a tiltrotor craft would yield a 35-percent power 
boost with a 6.5- to 7.3-percent improvement in Direct Operating Cost (DOC). The final design of the 
water- and methanol-injected engine had a smaller envelope and lighter weight with a 15-percent 
improvement in fuel bum for the 200 nautical mile mission over that of the baseline engine. 

Herein, we investigate combustor emissions for a trapped vortex combustor (TVC), as illustrated 
schematically in figure 1, and with diffuser and combustor hardware shown in figure 2. Liquid fuel is 
injected into the main combustion air through orifices in the individual tripass diffuser legs, and also into 
the upper and lower cavities. Water is injected directly into the upstream plenum, with sufficient 
residence time for complete evaporation. The plenum air and water flow feeds all combustion and cooling 
air injection ports and slots. Test rig instrumentation includes temperature and emissions probes along 
rakes at the combustor exit plane. Previous studies (Hendricks et al., 2001 and 2004) have characterized 
the TVC combustor flow fields, including detailed comparisons of predictions and measurements of 
aerodynamics and combustion performance. CFD code performance had been investigated and 
documented in a number of previous efforts (Ebrahimi et al., 2001; and Ryder and McDivitt, 2001) for 
different aerodynamic and reacting flow cases. In this effort, computed and measured results are 
compared for the case of cavity-only fueling. Specifically, comparisons are made of combustor exit 
temperature, CO and NO x emissions, at different water/fuel loading ratios. 


Trapped Vortex Combustor Rig 

The trapped vortex combustor (TVC) sector rig at the Air Force Research Laboratory, Wright- 
Patterson Air Force Base is used for the experiments. A schematic of the test rig is shown in figure 1, 
which illustrates the air and water flow and fuel injection sites, and the general flow patterns expected in 
the combustor. A photograph of the combustor hardware is shown in figure 2, with the side wall removed. 
Components of the rig include the tripass diffuser, combustor bulkhead, heat shield, and combustor duct 
that exhausts to a vent. Combustor walls are cooled through effusion holes along the entire interior wall, 
and film cooling through slots along the combustor upper and lower walls. Engine compressor exit air is 
emulated by a plenum, and is connected to the TVC combustor through a tripass diffuser (figs. 1 and 2). 
This diffuser configuration splits the flow evenly, with the upper and lower diffuser flow paths providing 
high-speed shear layers that drive the cavity flows and entrain products of combustion into the main flow. 
The center diffuser flow path provides a high-speed jet that interacts with the outer diffuser flows and 
further mixes the products as the flow approaches the combustor exit. The result is a stable combustion 
process with high efficiency, within a compact configuration. 

The TVC combustor is fueled in two different ways. Liquid fuel is injected through simple orifice 
injectors into the main diffuser flow paths, with the high-speed cross flow atomizing the liquid fuel jet 


NASA/TM— 2005-213441 


2 



into a fine mist. Independently, liquid fuel is injected into the TVC cavities through simple orifice 
injectors. The cavity spray droplets are injected into a hot, reacting flow environment and evaporate 
rapidly. The effects of the spray droplet injection conditions on the combustor flow field and combustor 
exit pattern factor has been investigated (Hendricks et al., 2004). During that effort, the initial droplet size 
distributions were determined for the engine operating conditions used in this investigation. For the main 
injector, and for a nominal plenum pressure of 50 psi (344.6 kPa), the Rosin-Rammler droplet size 
distribution was R 3 2 ,io = 3 pm, R 3 2 , 5 o = 6 pm, and R 3 2 , 9 o = 12 pm. For the cavity injector, the Rosin- 
Rammler droplet size distribution was R 32> io = 5 pm, R 39.30 = 8 pm, and R 32.90 = 10 pm. Through a 
combination of cavity-only, cavity and main, and main-only fueling schemes, different operating 
conditions can be produced in the rig, allowing for investigation of a wide variety of flow phenomena in 
the experiment. 


TVC Flow Visualization Studies 

The combustion flow in the TVC has been visualized in several ways, and examples are shown in 
figures 3, 4, and 5. Gray scale photos at different water loading ratios are shown in figure 3. With no 
water injection, the flame is highly luminous, with flame attached to the rear cavity wall. Entrainment of 
products into the main combustor stream is visible. At 50 percent water/fuel, the flame is less luminous, 
and the entrainment process is more difficult to observe directly. At 90 percent water/fuel, the flame 
appears faint, although combustion temperatures are only slightly lower. 

Time -averaged color photos for cavity-only fueling are shown in figure 4. On the left, combustion 
without water injection is shown, with a highly luminous flame filling the cavity, and entrainment into the 
combustor duct. A cool jet through the center is observed, due to the nonfueled combustor air entering 
through the central diffuser leg. On the right, flow with 90 percent water/fuel water injection is shown. 
The flame in this case is highly luminous, and appears better mixed than the nonwater-injected case. 

Still frames extracted from AVI video clips, taken at frame rates of 1 1 200 frames per second (fig. 5) 
show the variation in flame structure and luminosity over different conditions. In figure 5(a), a higher 
plenum pressure (100 psi or 689.3 kPa) is examined and displays a much higher luminosity than at lower 
pressures. Reduced fueling into the cavity (fig. 5(b)) results in a lower luminosity. Lower plenum pressure 
(50 psi, fig. 5(c)), and no water injection, still displays a high cavity luminosity, while 90 percent 
water/fuel water injection (fig. 5(d)) results in much lower luminosity. 

A detailed examination of the video clips show that at all pressures and water/fuel loading ratios, a 
dynamic combustion process, is taking place in the cavities. In particular, rotation of the primary vortex, 
precession of the vortices, and splitting of the primary vortex into two or more smaller vortices are 
observable. 


Flow Modeling 

Numerical simulations of the experiments were performed in order to determine whether the major 
trends in combustor exit quantities (temperature and emissions) could be predicted at the different water 
injection levels. To investigate this, a computational model of the TVC combustor rig was developed, 
including full geometry, and representation of all airflows. The model included the main diffuser airflow, 
cavity driving hole flow, cavity wall effusion cooling, slot cooling, and combustor wall liner film cooling 
(Hendricks et al., 2004). The combustion air and cooling airflow rates were determined through rigorous 
A cd analysis of the test rig, and were matched precisely in the simulations. 

In a related work (Hendricks et al., 2004), combustion was modeled using the eddy dissipation 
concept (EDC) heat release model, which proved adequate for predicting combustor exit temperatures, 
system pressure drop, and major flow structures with good agreement. In order to predict emissions in 
this effort, a hybrid chemistry model developed at the NASA Glenn Research Center, Cleveland, Ohio, by 
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Molnar and Marek (2003) was utilized. This hybrid model blends the mixing and ignition characteristics 
of the EDC model with the temporal accuracy of finite-rate chemistry. In developing the finite -rate 
chemistry formulation, hydrocarbon fuel kinetics equations for fuel and air combustion, including 
pressure dependency effects, were evaluated using stirred reactor codes. The computed results were used 
to develop kinetic reaction rates for a reduced three-step reaction mechanism: 


C\2H 2 3 +O 2 — 

XFuel >co+h 2 o 

(1) 

CO + j/,0 

2 ->CD 2 

(2) 

N 2 + O 2 - 

lNO >2NO 

(3) 


Chemical kinetic times for the fuel (x Fue i), CO (x C o), and NO x (x N o) are available as correlations 
covering a broad range of pressures, temperatures, and equivalence ratios, and are derived from GLSENS 
(see appendix) analysis of the 23-step, 16-species mechanism developed for Jet-A by Kundu et al. (1999), 
and Penko et al. (2000). Blending of the finite-rate component of the model with the EDC model is based 
on the evaluation of rate-limiting time steps. Computationally, switching between the EDC and the finite- 
rate model is dependent on the longer time scale determined for turbulent mixing rate and chemical 
kinetic formation rate (Marek and Molnar, 2003). 

Spray droplet initial size, droplet heat up and evaporation rate, and total droplet lifetime are key 
design parameters for injectors. These details must be understood for good matching of injectors and 
combustors to produce the desired heat release characteristics. In this effort, an analysis of the behavior of 
single droplets was performed to confirm aspects of the spray droplet modeling in the full CFD 
simulation. Specifically, single droplet characteristics were studied using a computer program developed 
for injection into crossflows (Liew et al., 2004). Injection into the cavities and diffuser legs were studied 
individually, with key results shown in figure 6. 

For the cavity injector, (figure 6, left column), spray droplets at the selected Rosin-Rammler sizes 
(5, 8, and 10 pm) and the computed liquid velocities were investigated for cold flow and hot flow ambient 
conditions. Cold flow ambient conditions represent a maximum droplet lifetime due to evaporation, while 
hot flow ambient conditions represent a minimum expected droplet lifetime. Inspection of the results for 
total droplet trajectory distance before complete evaporation (middle graph) indicate that the droplet 
trajectory analysis and the full CFD simulation results are in excellent agreement for the spray stream, 
with the stream fully evaporated before striking the opposite cavity wall. In particular, the smaller 
droplets evaporate more rapidly, while the larger droplets reside slightly longer in the cavity before 
reaching the boiling temperature and fully evaporating. 

For the main injector, the airflow is at a higher cross stream velocity; however, the analysis procedure 
is the same. Inspection of the right column of figure 6 shows similar trends as for the cavity injector, 
although the total distance traveled is much higher due to the high cross stream velocity. While large 
droplets injected into cold flow exit the combustor before complete evaporation, the hot flow condition 
shows these droplets fully evaporating within the combustor duct. This corresponds to the situation where 
some large droplets are not fully evaporated within the duct, and this is in fact the situation observed in 
the CFD simulations. 


Results 

Simulation cases were run to determine the overall trends in combustor exit plane temperatures and 
emissions for the case of main + cavity fueling. Following this, a set of cases was run corresponding to 
specific runs performed in the TVC rig at Wright-Patterson Air Force Base. For the numerical study, a 
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version of Run 2311 (Hendricks et al., 2004) was examined, utilizing liquid fueling into the mains and 
into the cavity. Water/fuel loadings studied included 0, 20, 40, 60, 80, 100, 150, and 200 percent 
water/fuel. Nominal plenum pressure was 50 psi (344.6 kPa). Water flow rates were additive to the 
airflow rates, keeping overall fuel/air ratios the same. Key graphical results from this numerical study are 
presented in figure 7, while mass-averaged quantities at the combustor exit are shown in figure 8. 

Computed results at the combustor midplane for fuel concentration, temperature, water, and CO mass 
fractions are shown in figure 7 at 20 and 200 percent water/fuel loading. In each graphic, instantaneous 
droplet positions are shown, with the droplet stream terminating at full evaporation. Cavity fuel droplets 
are seen to evaporate completely, while main fuel droplets streams strike the heat shield, reflect into the 
combustor duct, and generally evaporate before the combustor exit plane. Overall flow structure is similar 
for the 20- and 200-percent water loading cases. 

Mass-averaged combustor exit temperature initially decreases with water injection (fig. 8(a)), increase 
slightly at levels of 60 and 80 percent water injection, and then remain flat or decrease slightly at higher 
water loading ratios. Unbumed fuel at the combustor exit (sum of gaseous and liquid fuel), a measure of 
combustion efficiency, shows generally that as water loading is increased, more of the fuel is reacted. 
This may be partly due to the increased mixing observed in higher water loading cases. Carbon monoxide 
increases slightly as water loading is increased, while NO x decreases at the highest water loading, to a 
level of approximately one-half that observed without water injection. 

Graphical results for CFD predictions of a series of experimentally investigated cases at different 
water loading rates are shown in figure 9. Flow contours are shown at two water load levels 
(20 and 71 percent water/fuel); additional runs at 50 and 87 percent water/fuel were also performed. For 
the cases with cavity-only fueling, the cavity spray droplets are seen to evaporate rapidly, within the 
cavity. This changes little between the two water load ratios. The combustor duct flow is not as well 
mixed out in these cases, as in the cavity+main fueling runs (see fig. 7), due to the absence of fuel in the 
combustor center. 

A comparison of the post-processed, mass-averaged quantities at the combustor exit plane are shown 
in figure 10. Combustor exit temperature trends shown in figure 10(a) are similar for the predictions and 
measurements, with good agreement for the case without water injection. Increasing water fraction lowers 
exit temperature nearly linearly at levels up to about 75 percent water/fuel, with a steep drop of nearly 
200 K from the 71- to the 87-percent level. Predictions do not fully capture the magnitude of the 
temperature drop determined in the experiment, decreasing at a lower rate than in the measurements. CO 
levels at the combustor exit plane are shown in figure 10(b). Experiments and predictions both trend 
upwards with increasing water content, although the absolute prediction levels are significantly higher 
than the measured values. NO x levels at the combustor exit plane are shown in figure 10(c). Predictions 
and measurements both trend downwards, in a nearly linear manner, with the experiments showing a 
greater increase than predictions. 


Summary and Conclusions 

The effects of water addition on the combustion characteristics of the trapped vortex combustor 
(TVC) have been measured experimentally, and compared with simulations performed using a 
computational fluid dynamics (CFD) tool with a Jet-A chemistry model. Flow visualization indicated that 
flame luminosity decreased rapidly as a function of increasing water content in the combustion air. High- 
frame-rate videos of the flame revealed vortex dynamics in the TVC cavity, with vortex rotation, 
precession, and splitting into multiple vortex cores. For the investigated experimental runs, there was 
general agreement between the measured and computed trends in combustor exit temperature, CO, and 
NO x , although the model underpredicted the magnitude of the changes in each of these. The simulation 
results are being further investigated by using the collected experimental data, for the purpose of further 
evaluation of the Jet-A hybrid chemistry model in predicting trends in emissions levels for water-injected 
combustion processes. 
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Appendix 


GLSEN S is a Generalized Extension of LSENS (Bittker (1966)). LSENS is the NASA Lewis general 
chemical kinetics and sensitivity analysis code that computes the progress of complex multistep 
molecular processes in a homogenous gas mixture for several reaction models. For a static, nonflowing 
chemical reaction, it also enables the user to perform a sensitivity analysis of the reacting system 
simultaneously with the chemical kinetics computation. The sensitivity analysis is important to the CFD 
user, as it gives the user the ability to rapidly determine the relatively few individual reactions that are 
rate controlling in a system where many simultaneous steps are occurring for the static or quasi-static 
process (Radhakrishnan, Bittker (1994) in three parts and McBride and Gordon (1994, 1996)). 

Bittker, D.A.; GLSEN S, A Generalized Extension of LSENS, Including Global Reactions and Added 
Sensitivity Analysis for the Perfectly Stirred Reactor. NASA Reference Publication 1362, 1966. 
Bittker, D.A.; and Radhakrishnan, K.: LSENS- A General Chemical Kinetics and Sensitivity Analysis 
Code for Elomogeneous Gas-Phase Reactions, Part 111: Illustrative Test Problems. NASA RP-1330, 
1994. 

Gordon, S.; and McBride, B.J.: Computer Program for Calculation of Complex Chemical Equilibrium 
Compositions and Applications, Part 1-Analysis NASA RP-131 1, 1994. 

McBride, B.J.; and Gordon, S.: Computer Program for Calculation of Complex Chemical Equilibrium 
Compositions and Applications, Part 11-Users Manual and Program Description NASA RP-131 1, 
1996. 

Radhakrishnan, K.: LSENS- A General Chemical Kinetics and Sensitivity Analysis Code for 

Homogeneous Gas-Phase Reactions, Part 1: Theory and Numerical Solutions Procedures. NASA RP- 
1328, 1994. 

Radhakrishnan, K; and Bittker, D.A.: LSENS- A General Chemical Kinetics and Sensitivity Analysis 
Code for Homogeneous Gas-Phase Reactions, Part 11: Code Description and Usage. NASA RP-1329, 
1994. 
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Figure 1. — Schematic diagram of the trapped vortex combustor with water-injection strategy. Arrows 
indicating major flow components. Liquid fuel is injected into the TVC cavity and also directly into 
the main combustor air through orifices in the diffuser. Water mist mixed in with plenum air, feeding 
the cavity driving air and all cooling flows. 



Figure 2.— Test rig hardware, showing 10 fuel injector modules in the spanwise direction. 
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Figure 3.— Gray scale images of TVC, cavity fueling only, at different water-injection ratios, 
^cavity ~ 2-0 an d ^overall = 0-46. Left: 0 percent FLO. Middle: 50 percent water/fuel. Right: 
90 percent water/fuel. Plenum pressure = 50 psi (344.6 kPa). 



Figure 4.— Time-averaged color images of the flow. 4> Cavity = 2.2 and ()) 0verall = 0.46. Left: 

0 percent water/fuel. Right: 90 percent water/fuel. Both cases, plenum pressure = 50 psi 



100 psi, 4> Cavity = 0.8, and 0 percent water/fuel, (c) 50 psi (344.6 kPa), <t> Cavity = 2.2, and 0 percent 
water/fule. (d) 50 psi, (|) Cavity = 2.2, and 90 percent water/fuel. 
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Figure 6. — Single droplet trajectory and evaporation characteristics. Left column: cavity injector 
with R 32 = 5, 8, and 10 pm. Right column: main injector with R p = 3, 6, and 12 pm. Cold flow 
results using ambient injection conditions. Hot flow results obtained using averaged trajectory 
conditions. 
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Water and Fuel = 20% 


Water and Fuel = 200% 



Figure 7.— Cavity + main fueling. Midplane contours of combustion quantities, including evaporated 
fuel, static temperature, and water and carbon monoxide mass fractions. Spray droplets 
instantaneous positions shown. Left column: 20 percent water/fuel mass fraction. Right column: 
200 percent water/fuel mass fraction. 
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(a) Combustor exit temperature. 
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(b) Unburned fuel (gaseous + liquid) at 
mass fraction at combustor exit plane. 



(c) Carbon monoxide mass fraction. 
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(d) Nitric oxide (NO x ) mass fraction. 


Figure 8.— Mass-averaged quantities at combustor exit plane. Parametric study of water injection, 
ranging from 0 to 200 percent water/fuel mass fraction, with main and cavity injector fueling. 


NAS A/TM - 2005-21 344 1 


12 



Water and Fuel = 20% 


Water and Fuel = 71% 
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Figure 9.— Cavity-only fueling. Midplane contours of combustion quantities including evaporated 
fuel, static temperature, and water and carbon monoxide mass fractions. Spray droplets 
instantaneous positions shown. Left column: 20-percent water/fuel mass fraction. Right column: 
71 -percent water/fuel mass fraction. 


NAS A/TM - 2005-21 344 1 


13 


1800 


20 40 60 80 100 

Water/ Fuel Ratio (%) 


Combustor exit temperature. 
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(b) Carbon monoxide mass fraction at combustor exit plane. 
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(c) Nitric oxide (NO x ) mass fraction at combustor exit plane. 

Figure 10.— Mass-averaged quantities at combustor exit plane. Parametric study of water injection, 
ranging from 0 to 100 percent water/fuel mass fraction, with cavity-only injector fueling. 
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